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Fluorescence decay-time constants in organic liquid scintillators 
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The fluorescence decay-time constants have been measured for several scintillator mixtures based 
on phenyl-o-xylylethane (PXE) and linear alkylbenzene (LAB) solvents. The resulting values are 
of relevance for the physics performance of the proposed large-volume liquid scintillator detector 
LENA (Low Energy Neutrino Astronomy). In particular, the impact of the measured values to the 
search for proton decay via p K^V is evaluated in this work. 



I. INTRODUCTION 

Charged particles crossing a liquid-scintillator medium 
deposit their energy leading partly to fluorescence light. 
The lifetimes of the molecular excited states in the 
medium determine the pulse shape of the events in a 
detector application. For this reason, experiments re- 
lying on pulse shape analysis for background discrimina- 
tion are always dependent on the fluorescence decay-time 
constants. A common application of pulse shape analy- 
sis in neutrino experiments is the separation of a-events 
produced_by natural radioactivity from neutrino-induced 
/?-events [1] . The pulse shape also plays an important 
role in the search for proton decay via p K^V in the 
proposed large liquid-scintillator detector LENA as the 
signature relies on a double-peak structure in time . 
The LENA (Low Energy Neutrino Astronomy) detec- 
tor 3 is planned as a large- volume (50 kt) observatory 
based on the liquid-scintillator technology which will be 
highly suitable for the investigation of a variety of topics 
in astrophysics, geophysics and particle physics. In the 
analysis of the proton decay reaction p K^V, the two- 
peak time structure is used to discriminate proton decay 
events from the background due to atmospheric neutrino 
events Q . The detection efficiency depends upon the flu- 
orescence decay-time constants and therefore, the labo- 
ratory experiments presented here have been performed 
to determine these constants for several scintillator mix- 
tures. 

Fluorescence decay measurements of liquid scintillators 
have been performed earlier in feasibility studies for the 
Borexino experiment with a scintillator mixture based on 
PC (pseudocumene) as solvent In addition, similar 

measurements have been carried out for indium-loaded 
scintillators for the proposed solar neutrino experiment 
LENS [6] . A desirable scintillator for the proposed LENA 
detector would feature fast decay times, large attenuation 
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lengths and reasonable costs. This paper reports results 
obtained for the fluorescence decay constants of impor- 
tant current scintillator candidates. Non-metal loaded 
scintillators based on PXE and LAB have been tested. In 
addition, the impact of the decay constants on the sensi- 
tivity of the LENA detector to proton decay is shown. 



II. EXPERIMENTAL METHOD 

The main goal of the experiment is the determina- 
tion of the 'probability density function' (PDF) which 
describes the time-dependence of the photon emission in 
the fluorescence processes. A sampling of the photon- 
emission time after excitation of the medium is required. 
The experiment uses the 'start-stop' method [3] which 
measures the emission times for single photons. In this 
method, the start signal accounts for the starting point 
of the pulse: approximately the time of the energy de- 
position. If a light detector is able to measure a large 
number of photons of an event, the start time can be ex- 
tracted from the first detected photon. A second detector 
monitoring the same event provides the stop signal. This 
detector should be designed in such a way that the prob- 
ability of detecting a single photon of an event is in the 
few % regime. This guarantees that the probability to 
detect more than one photon in an event is almost zero. 
The PDF of a certain scintillator mixture results from the 
time difference between the start and the stop signals. 

The main setup consists of two photomultipliers 
(PMTs) which detect the light emitted from a liquid- 
scintillator sample. The sample is irradiated by a ^^Mn 
7-source (834 keV) . Figure[T] shows a scheme of the ex- 
perimental setup. A sample of scintillator is filled into a 
^ 10 ml cylindrical container made of black PTFE (Poly- 
tetrafluoroethylene, high density) material with quartz 
windows at both faces. The maximum energy for the re- 
coiling Compton electron for this source is 639 keV (back- 
scattering, i.e. Compton-cdge energy). Compton elec- 
trons deposit their energy in the scintillator leading to the 
fluorescence light-emission. One of the PMTs is directly 
coupled to the scintillator container (close PMT). It de- 
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tects a large number of photons per event and thereby 
sets the starting point for the event. The second PMT 
is placed ^ 30 cm away from the sample (far PMT) . This 
distance was chosen to assure single-photon events due 
to the small solid angle: The probability to detect a pho- 
ton from a Compton event in the far PMT is < 3%. This 
value has been calculated conservatively assuming a light 
emission of 10 000 photons/MeV 3 and an energy depo- 
sition corresponding to the Compton edge of ^''Mn. 
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FIG. 1: Schematic diagram of the setup for the fluorescence- 
time measurement. Two photomultiplier tubes (PMTs) col- 
lect the scintillation light emitted by a scintillator mixture. 
One PMT is directly coupled to the container and the second 
is located at a distance of 30 cm. The scintillator is excited 
by a ^''Mn 7-source. 



The whole setup is located inside a light-tight box. 
A NlM-based electronic setup has been used for signal 
processing. First, the signals of both photomultipliers are 
amplified and directly connected to the data acquisition 
system. The amplifier output for the close-PMT signal 
is, in addition, connected to a discriminator and used as 
trigger for the data acquisition system. 

For data readout and recording, an Acqiris (Ac- 
qiris/ Agilent Technologies Q ) data acquisition system is 
used. In the setup, the acquisition system is used at 
a sampling rate of 4GS/s (giga-samples per second) to- 
gether with a 10 bit dynamic range. The system is op- 
erated in 'external trigger mode' with the close PMT 
providing the hardware trigger-signal. A Labview-based 
software is used to generate a software trigger and to 
store the data. When a trigger signal from the close mul- 
tiplier appears, a coincidence window of 500 ns is opened. 



If the software detects a signal above threshold in the far 
PMT, both photomultiplier data pulses are stored. The 
resulting files contain 50 000 pulses covering a measuring 
time of 8h. 

The time resolution of the system has been measured 
by a slight modification of the setup just described. The 
^^Mn source is replaced by a ^^Na /3+ source which pro- 
duces two 511 keV 7-rays from e"*" e~ annihilation. The 
far PMT is optically decoupled from the fluorescence sys- 
tem by a black blocking filter. A Cherenkov radiator 
(Plexiglas piece) is placed in front of the far PMT. This 
provides an instantaneous signal as the Cherenkov pho- 
tons are produced within picoseconds time after the en- 
ergy deposition of a 511 keV gamma. Pairs of 7-rays 
from the annihilation are emitted isotropically but al- 
ways in opposite directions. Sometimes both gammas 
interact simultaneously in the liquid scintillator and in 
the Cherenkov radiator and those coincident events have 
been recorded. Using the time difference between one 
511 keV signal in the Cherenkov radiator and a second 
one in the liquid scintillator sample, the time resolution 
of the system is determined. 



III. TIME RESOLUTION AND BACKGROUND 
SOURCES 

The time resolution of the fluorescence-measurement 
system mainly depends upon the single-photoelectron 
transit time jitter of the PMT. Photoelectrons produced 
at different positions of the cathode have different propa- 
gation distances to the first dynode resulting in a smear- 
ing of the transit time through the tube. The specifica- 
tion for the time jitter of the PMTs used is 0.98ns[9|. 
Smaller contributions to the loss in time-resolution are 
the error of the determination of the start time of a pulse 
and the time jitter of the electronic readout chain, the 
latter, however, usually being a small contribution. 

FigurelH shows the distribution of the time differences 
between the start signal given by the scintillator and the 
arrival of the Cherenkov photons measured in the setup 
with the ^^Na source. Although the main peak is well 
described by a Gaussian function, there are small contri- 
butions at times below and above the main peak. The 
shoulder on the left side of the Gaussian curve is due to 
pre-pulses of the single-photon tube. Further to the left, 
there is a contribution which can be interpreted as pho- 
tons which cross the cathode producing a photoelectron 
directly in the dynode-chain. The total contribution of 
the events to the left side of the Gaussian peak is about 
5% of all pulses. 

The pulses on the right side of the Gaussian peak 
(t < 20 ns) can be explained by the reflection of photo- 
electrons by the first dynode back to the cathode. This 
effect occursjn about 2% of all pulses and has also been 
reported in 4| . The pulses at still later times {t > 20 ns) 
contribute about 2% to the total number of pulses. Their 
origin is unknown but may be interpreted as a weak flu- 
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FIG. 2: Distribution of the time differences between the start 
signal in the scintillator and the arrival of Cherenkov photons. 
The time resolution of the system (the events are plotted on a 
semi-logarithmic scale) is well described by a Gaussian curve. 



orescence of the photomultiplier housing and window as 
the radioactive source was placed very close to the tube. 
When removing the Plexiglas piece, it was found that the 
number of pulses with f > 20 ns increase relative to the 
previous measurement pointing to a background related 
to direct 7-hits on the PMT. Altogether, about 91% of 
the pulses are beneath the Gaussian curve. 

In addition to the loss of time resolution due to the 
effects described above random coincidences due to the 
dark-count rate of the far PMT contribute to the back- 
ground of the experiment. This coincidence rate has been 
estimated by measuring the dark current of the far PMT 
and taking the trigger rate by the close PMT. About 
30 random coincidences are expected for an 8 hours mea- 
surement 0.06% of all recorded pulses). 



IV. MEASURED SAMPLES 

The samples investigated consist of a mixture of an 
organic liquid solvent and one or two wavelength-shifters 
as solutes. The solvents PXE (phenyl-o-xylylethane) and 
LAB (linear alkylbenzene) have been studied because of 
their high light yield and low risk properties. 

The wavelength-shifters dissolved were PPO (2,5- 
diphenyl-oxazole) , pTP (para-terphenyl), bisMSB (1,4- 
bis-(o-methyl-styryl)-benzene) and PMP (l-phenyl-3- 
-mesityl-2-pyrazoline). To characterize these compounds 
a spectrum of the emitted light has been recorded when 
the compounds were dissolved in PXE. To this end, 
the medium was excited with UV-light of a deuterium 
lamp [lO| . The spectra were recorded with a compact 
spectrometer with 0.53 nm resolution and fiber optic cou- 
pling (Ocean Optics IIR2000) [ll|. They were corrected 
with the spectrometer's sensitivity curve and also for pos- 
sible shifts in the wavelength |12| . Figure[3] shows the 
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FIG. 3: Emission spectra of PPO, bisMSB and PMP, all with 
a concentration of 2 g/£ in PXE solvent. 



emission spectra of the wavelength-shifters PPO, bisMSB 
and PMP. The data is normalized such that the max- 
imum of the spectrum is at 1. All solutes have been 
dissolved with a concentration of 2 g/£ in PXE. 
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V. ANALYSIS AND RESULTS 

The first step in the analysis is the determination of 
the start time for both recorded pulses (close and far 
PMTs). The start time is derived in a constant- fraction 
method by setting the starting point of each pulse to the 
40% height in the leading edge of the pulse. The distribu- 
tion of the photon emission time results from the differ- 
ence between the starting point of the close PMT pulse 
(zero-point) and that of the far PMT. Pulse height cuts 
are performed to reject dark-current pulses and overflow 
pulses below and above the single photon peak, respec- 
tively. In addition, pulses with more than one peak in 
the single-photon PMT are also rejected (double pulses 
or noise pulses). 

The probability density function (PDF), F{t), of the 
scintillation-light emission can be described by a convo- 
lution of a multi-exponential function with the time res- 
olution of the measuring system: 




where is the decay-time constant of the exponential 
decay function i and Ni is the contribution of this com- 
ponent such that Ni = 1; to is the time at which the 
exponential decay starts and R{t) is the time-resolution 
function. In good approximation, the time resolution can 
be represented by a Gaussian curve: 

R{t) ^ (2) 

where cr is a known parameter. The time resolution of 
the system is affected by the uncertainty in the deter- 
mination of the pulse start-time mentioned above. For 
the close PMT, the start time is more accurate for scin- 
tillators with a high light yield due to the better pho- 
ton statistics. Therefore, the resolution depends on the 
scintillator sample and has to be determined for every 
mixture. This has been achieved by fitting a Gaussian 
curve to the left (sharply rising) part of the measured 
PDF distribution. This method gives compatible results 
with the direct time-resolution measurement described in 
section lllll 

The ROOT-based Roofit toolkit is used to per- 
form the fit of equation[T] to the measured data. In 
its default configuration, Roofit performs a maximum 
likelihood fit and minimizes — log(L) where L is the 
likelihood function. This minimization is performed by 
calling the Minuit package from CERN libraries. 
FigureH] shows the measured PDFs for four examples 
of different scintillator mixtures: PXE -I- 2g/^PP0 -I- 
20mg/£bisMSB, PXE + 2g/£PP0, LAB + 2g/£PP0 
and PXE + 2g/^PMP. The data is presented in his- 
tograms with a bin size of 1 ns together with the statis- 
tical errors. A fit using the convolution of the sum of 
several exponential functions with a Gaussian curve is 



TABLE I: Nomenclature used for the mixtures of which the 
fluorescence decay constants have been measured. The sub- 
index c corresponds to the concentration of the corresponding 
solute. 



Name 


Mixture 


PPc 


PXE + cg/£PPO 


LPc 


LAB + cg/iFPO 


PMc 


PXE -1- cg/£PMP 


PPb 


PXE+2 g/^ PPO+20 mg/e bisMSB 


PTb 


PXE+2 g/£ pTP+20 mg/i bisMSB 



also depicted by the red curve. 

The abbreviated names of the mixtures for which the 
fiuorescence decay-time constants have been measured 
are summarized in table|Tl The sub-index c corresponds 
to the concentration of the corresponding solute. 

TablesHI] - |V] summarize the parameters of the fits ac- 
cording to equation[T] Depending on the scintillator mix- 
ture and the solute concentration, either three or four 
decay components perform the best fit. Fit parameters 
are the time resolution cr, the decay constants r^, the 
amplitudes Ni of each exponential function and the start 
time to of the decay function. The amplitude of the last 
component {Ni^cah third of fourth depending on the fit) 
has been calculated according to J^i — ^ using the 
fit values. Its error depends on the other amplitudes. 
To estimate a maximum error [ANmax) for this param- 
eter, quadratic error propagation has been used. The 
uncertainties Ar^ and ANi consist of a statistical error 
Sstat which is obtained from the fit and a systematic er- 
ror, Ssys- The latter has been derived from a systematic 
study where the parameter a is varied within one stan- 
dard deviation of its value (a + Aa, a — Aa). It was found 
that the systematic error is of the same order of magni- 
tude as the statistic one. Both systematic and statistic 
uncertainties are in general asymmetric. To evaluate the 
quality of a certain fit, both the convergence of the fit 
and the value of the likelihood function (— log(_L)) have 
been considered [10|] . 

VI. DISCUSSION 

For each sample, the shortest decay-time constant ob- 
tained, Ti, has a value of few nanoseconds (2 — 8 ns) and 
represents the main contribution of the photon emission 
(its amplitude varies between 60% and 95% of the total 
emission). This contribution is related [l^ to the transi- 
tion between the lowest excited singlet spin-state 5*1 in 
the solute and its ground state ^o. The additional two 
or three components arise from further molecular pro- 
cesses such as the de-excitation of electrons in triplet 
spin-states. In general, the data shows that mixtures 
based on LAB solvent show longer decay constants. This 
observation has been made by other groups [l3| as 
well. Concerning the wavelength-shifters, PPO shows 
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FIG. 4: Probability density function of the scintillator mixtures (mix) PXE + 2g/^ PFO + 20 mg/^ bisMSB, FXE + 2g/.^FFO, 
LAB + 2g/^PP0 and FXE + 2g/i?FMF. The plots show the time dependence of the light emission fitted (curve in red) 
according to equation[T] 



TABLE 11: Resuhs of the fits for the mixtures based on PXE + 2 g/^ PPO and 2 g/^ pTP, both with 20 mg/^ bisMSB as 
secondary wavelength-shifter. All decay times Ti, the time to, as well as the time resolution a are given in ns. The amplitudes 
Ni are given in % with Ni^cai calculated from Ni = 1. The errors Ar; and AA''i include both statistic and systematic 
uncertainties. 



Mixture (mix) 


PPb 


PTb 


a (ns) 


1.61 ±0.01 


1.50 ±0.01 



mix 


Ti ± Ari 


T2 ± At2 


T3 ± At3 


Ti ± Ar4 


to ± Ato 


PPb 


^•"-■^-0.02-0.03 


7 7+U.2+U.2 
' • ' -0.3-0.3 


0/I+2+U 
•J^-2-1 


205ti^ti 


a 7rq+0.UU6+U.01 

U. (dO_Q QQg_Q Ql 


PTb 


r, qo+0.04+0.05 
^••J°-0.04-0.04 


7 0+0.2+0.1 
' •'J-0. 2-0.1 


oc+:i + i 

•J^-2-O 




0-00^-0.008-0.02 



mix 


iVi ± ANi 


N2 ± AiVa 


Ns ± AiVg 




PPb 


n 8c;q+0-007+U.005 

U.OdO_Q QQg_Q QQg 


n nQ7+**-U07+0.0U5 
u.uy r -0.006-0.005 


U.UO^_0. 002-0 


0.02 ±0.01 


PTb 


r, f.r,+U.01+0.U3 

'-'•O^-O.Ol-O.Ol 


n qn+D.Ol+O 
O-'J'J-O.Ol-O.Ol 


n r,C7+U.U02+U.003 
U.UO/_o. 002-0 


0.02 ±0.01 



shorter values for ti than PMP. The addition of the sec- 
ondary wavelength-shifter bisMSB to a mixture contain- 
ing PPO has only a small effect on the decay times. 

For three scintillator systems (PXE-PPO, LAB-PPO 
and PXE-PMP), the evolution of the shortest decay time 



constant ti with the solute concentration c has been stud- 
ied. From the variation of this decay-time constant, in- 
formation on the energy transfer in the different systems 
can be obtained. The energy transfer mechanism can be 
described in the following way [l^ Q . After the energy 
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TABLE III: Results of the fits for mixtures of tlie solvent PXE with different concentrations of PPO. Further details as in 
tableim 



mix 


PPl 


PP2 


PP6 


PPIO 


a (ns) 


1.45 ±0.01 


1.229 ±0.007 


1.066 ±0.006 


1.171 ±0.007 



mix 


n ± An 


T2 ± Ar2 


T3 ± At3 


r4 ± Ar4 


to ± A<o 


PPl 


o -|f.+U.IM+U.U6 
"-••^"-0. 05-0. 06 


7 7+o.v+o.y 

' ■ ' -0.7-0.5 


34+^+^2 


910+^5+12 


a (^01 +0.006+0.01 

U.U01_Q QQg_Q QJ 


PP2 


9 (.0+0.02+0. 01 
^•"'-'-0.02-0.02 


o O+0.3+0.1 
"•"-0.3-0.2 


40+3+0 


0/1 0+*^ +"^ 
^^^-21-1 


a Rf.-, +0.005+0.007 
"•""^-0.005-0.008 


PP6 


9 r)Q+0.01+0.01 

^•'^'3-0.01-0.01 


q Q+0.3+0.1 

•^•'^-0.3-0.1 


^7+3+0 
^' -3-1 


9r,o + 16 + l 

^"•3-13-0 


fi 7Q'>+0. (104+0.006 
"•' '3 '3-0. 004-0. 007 


PPIO 


9 n/l +0.01+0.02 
^•'-'^-0.01-0.01 


o f.+0.4+0.2 
°-"-0. 4-0.1 


40i^i^3 


207+14 + 1 


a R97+0. 005+0. 008 
U.UZ I -0.005-0.008 



mix 


iVi ± ATVi 


N2 ± A7V2 


iV3 ± ATVg 




PPl 


n +0.02+0. 03 

"•°^-0.02-0.02 


n -,9+0.02+0.02 

"•-■^^-0.01-0.02 


"•"^^-0.002-0.002 


0.011 ±0.028 


PP2 


n 0C9+0. 005+0. 002 

"•""-"^-0. 005-0. 003 


n in'?+0-005+o.oo2 

"•-■^""3-0.004-0.003 


n n-?! +0.001+0 

"•""3^^-0.001-0 


0.014 ±0.007 


PP6 


r, 099 + 0.003+0.001 

"•'5^^-0.003-0.002 


n -inc;+o. 002+0. 001 

"•-•^"^-0.002-0.001 


"•"^"-0.002-0 


0.027 ±0.004 


PPIO 


n ono+o. 004+0. 002 

"•""^-0.004-0.002 


"•-•""3-0. 003-0. 002 


"•""^-0.002-0.001 


0.032 ±0.005 



TABLE IV: Results of the fits for mixtures of the solvent LAB with different concentrations of PPO. Further details as in 
tableHI] 



mix 


LPl 


LP2 


LP6 


LPIO 


a (ns) 


1.92 ± 0.01 


1.58 ±0.01 


1.236 ±0.009 


1.133 ±0.008 



mix 


n ± An 


T2 ± At2 


T3 ± At3 


r4 ± Ar4 


to ± Ato 


LPl 


7 4g+o +0.04 

' •^"-0. 07-0. 03 


99 0+0.6+0.2 
^^•"3-0.6-0.2 


115t^±^i' 




A p,-, +0.01+0.01 
^•"-•--0.01-0.01 


LP2 


r 9-1+0.04+0.03 
"'■^-'--0.05-0.03 


18 /I+0.6+0.2 
-'^^•4-0.6-0.2 


118+^+'i^ 




A ifii +0.009+0.009 
^■^"J^-0. 009-0. 01 


LP6 


9 7-1+0.06+0.07 
^- ' -"--0.08-0.08 


7+0.6+0.5 
"• ' -0.6-0.4 


30±^±i 


136+^i 


A C77+0.007+0.01 
^■"^' ' -0.007-0.01 


LPIO 


1 Q/1 +0.04+0.04 
-^•■^^-0.05-0.04 


r Q+0.3+0.2 

"-"••3-0.3-0.2 


^"••^-0.9-0.3 




A ccn+O.UOV+O.Ol 
^•"^"3"-0. 007-0. 01 



mix 


iVi ± AiVi 


N2 ± AiV2 


iV3 ± AiV3 




LPl 


n 7c;q+(). 009+0. 004 
"• '"3 -3-0. 009-0. 003 


n 9-1 +0.008+0.003 
"•^-•--0.008-0.003 




0.031 ±0.012 


LP2 


n 777+0.007+0.003 
"• ' ' ' -0.007-0.003 


n 1 7n+('-"0t'+0-tJ(>^ 

"•3 ' "-0. 006-0. 003 




0.053 ±0.009 


LP6 


n f;c+0.03+0.03 
"■""3-0.04-0.03 


n 9-, +0.02+0.03 
"•^^-0 -0.02 


n inn+"-U04+o.oo2 

"■3""-0. 004-0. 003 


0.040 ±0.036 


LPIO 


n r,K+0.02+0.01 
"•"3"-0.02-0.01 


r, 97+0.01+0.01 
"•^ ' -0.01-0.01 


n -, 00+0.004+0.01 
"•^JO-o. 004-0. 02 


0.037 ±0.023 



TABLE V: Results of the fits for mixtures of the solvent PXE with different concentrations of PMP. Further details as in 
tableHIl 



mix 


PMl 


PM2 


PM6 


PMIO 


a (ns) 


2.12 ±0.01 


1.77 ±0.01 


1.362 ±0.009 


1.258 ±0.009 



mix 


n ± An 


T2 ± At2 


T3 ± At3 


to ± Ato 


PMl 


A qn+O.Ol+O.Ol 
^■"3"-0.01-0.02 


9-1 .+0.6+0.2 
^3-'i_o.6-0.2 


91 9+22+3 


A /ifiR+0. 007+0. 009 
^•^""-0.007-0.009 


PM2 
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deposition by a charged particle, the solvent molecule 
(PXE or LAB) rapidly reaches the first excited state. 
Via collisions with neighboring solvent molecules, the ex- 
citation energy can propagate between solvent molecules 
with a rate kt [s~^]. This process is also known as energy 
hopping [l9l |. The excitation energy moves spatially un- 
til a solute molecule is encountered to which the energy 
is transferred by dipole-dipole interaction 20] . The mea- 
sured decay constant ti is mainly the sum of the energy- 
hopping time Tt (c) = ^^^^ and the intrinsic lifetime of the 
solute Ts, Ti = Ts + Tt{c). The number of solvent-solvent 
collision processes n depends on the solute concentration 
c, the higher the concentration the lower the value for n, 
the faster the energy transfer from a solvent to a solute 
molecule. 

Assuming a linear dependence of the transfer rate on 
the solute concentration [21] [22] and defining an effective 
hopping rate kh, it follows that: kt/n{c) — kh ■ c/c^ and 
the measured decay constant ri can be written as: 



n c 



'ih- c' 



(3) 



where cq — 1 g/^. By fitting this formula to the exper- 
imental data, the intrinsic lifetime of the solute can 
be extracted. Moreover, the value of k^ represents the 
efficiency of the energy transfer via collisions (hopping) 
between solvent molecules until the energy reaches a so- 
lute molecule. 

Figure[5] shows the evolution of ti with the PPO con- 
centration for PPO being dissolved in the solvents PXE 
(upper panel) and LAB (middle panel) . The lowest panel 
shows the concentration dependence for the solute PMP 
dissolved in PXE. The errors assigned to the data points 
are the linear sum of the statistic and systematic uncer- 
tainties. The fit applied to the data points corresponds 
to equation[31 For each plot, the right upper-box shows 
the fit results where the parameter is labeled with the 
name of the solute and kh with the name of the solvent. 
The fit takes both the uncertainties in the solute concen- 
trations and the errors in the decay-time constants into 
account, with the latter errors being asymmetric. 

From the fits, the intrinsic lifetime of the PPO singlet- 
state Si can be extracted. For the PXE-PPO mix- 



tures T, 



= 1.85 ± 0.03 ns is obtained and r. 



1.19 ± 0.15 ns for the mixtures of LAB-PPO. Thus, the 
obtained value for Ts is lower when the LAB solvent is 
used. The solute lifetime changes with the solvent type 
because the interaction between both molecules changes. 
In the literature, the value for the intrinsic lifetime of 
the Si state of PPO varies slightly. In [2^ the values 
1.44 and 1.8 ns are reported. They result from different 
theoretical calculations of the Si Sq transition and in- 
clude a factor which takes the index of refraction of the 
solvent medium into account (cyclohexane in the pub- 
lished cases). Finally, in the value 1.6 ns has been 
used. The measured values for the PPO lifetime do not 
agree within the given error. We attribute this differ- 
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FIG. 5: Evolution of the fluorescence decay time ri for differ- 
ent concentrations of PPO in the solvents PXE (upper part) 
and LAB (middle). Bottom: evolution of ti for different con- 
centrations of PMP in PXE. 



ence to solvent effects. The theoretical values given in 
the literature are in fair agreement with the results for 
Ts presented here. 

The parameter kh quantifies, as explained above, the 
transfer of energy between solvent molecules via colli- 
sions. This transfer rate is faster for PXE than for LAB 
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which explains why for low PPO concentrations the mea- 
sured Ti values are much shorter than those for PPO 
dissolved in LAB. As an example, for a concentration 
of 2g/^ PPO, Tt_pxE = 2.63ns and Tt.LAB — 4.03ns (see 
equation[3]). 

The obtained value for the intrinsic lifetime of PMP, 
Ts.PMP = 3.30 ± 0.05 ns is close to the measured value 
of 3.01ns given in reference . Ts^pmp is significantly 
larger than the intrinsic lifetime of PPO. This might 
be due to the fact that the overlap between absorption 
and emission spectra in PMP is significantly lower than 
in PPO. The value of fc/i,pxE of PMP is compatible to 
the one measured in the PXE-PPO system (fc^^pxE — 
0.57 ± 0.09 when PMP is used and fcft,pxE = 0.72'± 0.08 
for the mixture with PPO). The measured ti constants 
decrease with increasing concentration and can well be 
fitted by equation[3l However, for large detectors solute 
self-absorption increases with concentration as well lead- 
ing to longer ti . This is due to the additional time which 
the light reemission takes [l^l ■ 

VII. IMPACT ON SCINTILLATION-DETECTOR 
PERFORMANCE 

The measurements of the fluorescence decay-time con- 
stants show that the main component of the light emis- 
sion (ti , A^i ) is significantly dependent on the scintillator 
mixture. The impact of this variation on the potential 
of the proposed LENA detector to search for proton de- 
cay {p — > AT+F) has been investigated in the present pa- 
per. The signature of these decay events consists of two 
pulses close in time f2f . The first arises from the energy 
deposited by the kaon, the second from the energy of 
its decay particles (Ar+ ^ /i+ + and — > 7r+ -I- 7r°, 
with 63.4 and 20.9% branching ratios, respectively). Due 
to the short lifetime of the kaon (12.4 ns) these signals 
overlap significantly. However, the double-peak struc- 
ture in the time can be used to discriminate those events 
from atmospheric neutrino background. The sensitivity 
of LENA to search for proton decay has been studied with 
a Geant4 l2y|-based Monte Carlo simulation presented in 
reference [2[ . In those simulations, a main decay time 
constant ri of 3.5 ns was used. A sensitivity to the proton 
lifetime of 4 • 10'^'* y has been derived from the analysis. 

Using the results of the present work, new simulations 
have been performed for two possible values of n (3 and 
6 ns) which are reasonable according to the presented 
measurements. For some samples, the decay time con- 
stant Ti measured was even shorter than the chosen val- 
ues. However, light propagation also has an influence on 
the pulse shape as scattering and solute self-absorption 
smear slightly the time information. The values of the 
absorption and scattering lengths are set to 20 m each 
as laboratory measurements have shown that those are 
realistic values. The resulting attenuation length is 10 m. 
All other scintillation constants were kept at their previ- 



ous values. The time jitter of the photomultipliers is also 
included in the simulation. Further details are discussed 
in 0. The resulting efficiencies for proton decay detec- 
tion in LENA via p AT+F are 65% and 56% for the 
decay values ti = 3 ns and ri — 6 ns, respectively. This 
means that the choice of the scintillator mixture can have 
an effect of about 20% on the achievable proton lifetime 
limit. 



VIII. CONCLUSIONS 

Fluorescence decay-time constants have been studied 
for several scintillator mixtures which could be relevant 
for the LENA detector. The effects of variation of the 
solvent type, solute type and concentration have been 
investigated systematically. In the measured samples, 
there is a short decay-time constant (ti) which accounts 
for the main part of the emitted light (60 to 95%). For 
this decay constant, large differences from 2 to 8 ns have 
been obtained for various scintillator mixtures. Concern- 
ing the fluorescence decay-time constants, attenuation 
lengths and the solute self-absorption, a mixture con- 
taining PXE, about 2 g/e PPO and - 20 mg/i bisMSB 
shows the best performance for LENA. The emission 
spectrum of bisMSB is shifted to longer wavelength val- 
ues (see figure|31) in comparison with spectra of usual pri- 
mary wavelength-shifters like PPO. Since absorption and 
scattering lengths increase with the wavelength, the ad- 
dition of bisMSB is favored. In principle, mixtures of 
PXE with PPO show a shorter ti but the self absorption 
of the wavelength-shifter plays also a role in the light 
propagation over large distances smearing the time in- 
formation. The addition of a mineral oil like dodecane is 
currently under discussion. It has the advantage of be- 
ing very transparent (long attenuation length) and has a 
large number of free protons. The latter is important for 
the detection of via their capture on free protons [l^ . 
However, the addition of dodecane would result in slower 
decay-time constants [2^ . 

The performance of the liquid scintillator in LENA 
has direct implications on the physics potential. This 
has been studied within a Monte Carlo simulation. If 
the shortest fluorescence decay-time constant is increased 
from 3 to 6 ns, the detection efficiency and therefore the 
sensitivity is reduced by 20%. 
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